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Parasite load and risk factors for poor outcome among 
children with visceral leishmaniasis. A cohort study 
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Clinical and laboratory risk factors for death from visceral leishmaniasis (VL) are relatively known, but quan- 
titative real-time polymerase chain reaction (qPCR) might assess the role of parasite load in determining clinical 
outcome. The aim of this study was to identify risk factors, including parasite load in peripheral blood, for VL poor 
outcome among children. This prospective cohort study evaluated children aged < 12 years old with VL diagnosis at 
three times: pre-treatment (TO), during treatment (TT) and post-treatment (T2). Forty-eight patients were included 
and 16 (33.3%) met the criteria for poor outcome. Age < 12 months [relative risk (RR) 3.51; 95% confidence interval 
(CI) 1.89-6.52], tachydyspnoea (RR 3.46; 95% CI 2.19-5.47), bacterial infection (RR 3.08; 95% CI 1.27-7.48), liver 
enlargement (RR 3.00; 95% CI 1.44-6.23) and low serum albumin (RR 7.00; 95% CI 1.80-27.24) were identified as 
risk factors. qPCR was positive in all patients at TO and the parasite DMA was undetectable in 76.1% of them at Tl 
and in 90. 7% at T2. There was no statistical association between parasite load at TO and poor outcome. 
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It is estimated that there are between 200,000- 
400,000 cases of visceral leishmaniasis (VL) in the 
world annually. Over 90% of these occur in six coun- 
tries: India, Bangladesh, Sudan, South Sudan, Ethiopia 
and Brazil (Alvar et al. 2012). In Brazil, there has been 
an increase in the number of the VL cases over the past 
few decades, with an annual average of 3,484 cases be- 
tween 2000-2010 (MS 2011a, b). Over the last decade, a 
high lethality rate has been observed, ranging from 3.2- 
8.5% (MS 2011a, b). 

Belo Horizonte, the capital of the state of Minas 
Gerais, an endemic region for VL, also showed an in- 
crease in VL cases, with an annual average of 83 cases 
per year between 2000-2004 and 132 cases per year be- 
tween 2005-2010. The lethality rate in Belo Horizonte 
between 2006-2011 was higher than the national aver- 
age, ranging from 8.1-23% ( PBH/SMS 2012). 

Previously identified risk factors for death are pro- 
longed course of the disease, very old and very young 
age, malnutrition, vomiting, diarrhoea, jaundice, bleed- 
ing, dyspnoea, oedema, coinfection with human immu- 
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nodeficiency virus, bacterial coinfection, severe anaemia, 
leucocytosis, neutrophilia and neutropaenia (Seaman et 
al. 1996, Lyons et al. 2003, Werneck et al. 2003, Collin et 
al. 2004, Rey et al. 2005, Mueller et al. 2009, Costa et al. 
2010, Sampaio etal. 2010). 

Various infectious diseases, such as acquired immune 
deficiency syndrome and hepatitis C, have a clear corre- 
lation between microorganism load and disease outcome 
(Ghany et al. 2011, Panel on Antiretroviral Guidelines 
for Adults and Adolescents 2013). Two studies assessed 
the association between pre-treatment parasite counts 
and VL outcome. These studies, from Sudan and Brazil, 
showed that higher number of parasites in the spleen or 
bone marrow (BM) were associated with a higher death 
rate (Seaman et al. 1996, Costa et al. 2010). Both studies 
used microscopic methods for parasite quantification, 
which have limitations due to the invasive procedures 
required to obtain the spleen and BM aspirates and the 
excessive time required for the analysis of these samples 
(da Silva et al. 2005). Quantitative real-time polymerase 
chain reaction (qPCR) is a technique that can quantify 
parasite load in an accurate manner in different biologi- 
cal samples, including peripheral blood. Moreover, it pro- 
vides results more rapidly than microscopic techniques, 
with less DNA contamination and fewer false positive 
results than standard PCR (Kubista et al. 2006). 

Considering VL severity and lethality rate, the early 
identification of factors associated with poor outcome 
and death is relevant to the premature establishment of 
appropriate treatment. The objective of this study was to 
identify risk factors, including parasite load quantified 
by qPCR in the peripheral blood, for VL poor outcome 
among children. 
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SUBJECTS, MATERIALS AND METHODS 

Patients and clinical assessment - This was a pro- 
spective cohort study conducted at Joao Paulo II Hospi- 
tal (HJPII) of Hospital Foundation of Minas Gerais State 
(FHEMIG) from June 2010-June 2011. 

The inclusion criteria were hospitalisation at HJPII, 
age < 12 years old and laboratory-confirmed VL by one 
of the following methods: BM smear, PCR (k-DNA) in 
BM aspirate or in peripheral blood, indirect immunofluo- 
rescence antibody test or rk39 immunochromatographic 
antigen test (Kala-azar Detect™ or Diamed IT-Leish™). 
All patients with confirmed VL in the study period were 
invited to participate the study. 

The main endpoint was poor outcome, characterised 
by at least one of the following interventions: admission 
to the intensive care unit (ICU), administration of vaso- 
active amines, mechanical ventilation or transfusion of 
blood products (red blood cells, platelets or fresh frozen 
plasma). Demographic and clinical characteristics and 
laboratory findings were assessed before the initiation 
of treatment. 

The patients underwent three clinical evaluations: 
pre -treatment (TO), between 10-15 days of treatment 
(Tl) and between 40-60 days after initiation of treatment 
(T2). Spleen and liver sizes were measured, considering 
their extensions below the costal margin at midclavicu- 
lar line. Then, the measures were divided by child body 
surface, in order to obtain an adjusted measure, compa- 
rable between all children, despite the weight and height 
variation observed among these patients. In each clinical 
evaluation, peripheral blood samples were withdrawn 
for parasite load assessment. For the continuous mea- 
sures, such as laboratory values, cut points were defined 
according to Brazilian and international guidelines. 

Parasite load by qPCR quantification - Total DNA 
from peripheral blood samples was extracted using the 
QIAamp DNA Mini Kit (Qiagen GmbH, Hilden, Ger- 
many), according to the manufacturer's instructions. Two 
independent assays for the detection and quantification of 
Leishmania spp and human DNA were performed using 
the StepOnePlus™ Real-Time PCR System (Life Technol- 
ogies, Carlsbad, CA, USA). For the Leishmania assay, the 
target DNA was the small-subunit ribosomal RNA (SSU 
rRNA) gene, which is conserved among all Leishmania 
species. The Leishmania assay consisted of the primers 
LEIS.U1 (5'-AAGTGCTTTCCCATCGCAACT-3') and 
LEIS.L1 (5'-GACGCACTAAACCCCCTCCAA-3'), de- 
signed to amplify a 67-bp fragment and the fluorogenic 
probeLEIS.Pl (FAM5'-CGGTTCGGTGTGTGGCGCC- 
3'TAMRA). The protocol described by Gomes et al. 
(2012) was applied. For the human assay, the ACTB refer- 
ence gene was used as the target and the primers Acol 
and Aco2 (Musso et al. 1996), which generate 120-bp 
fragments, were used in this assay. The reaction mixtures 
contained 12.5 uL of Syber® Green PCR Master Mix 2X 
(Life Technologies), 0.1 jj.M of each primer and 3 uL of 
DNA template in a final volume of 25 uL. The cycling 
parameters were universal and the melting analysis was 
conducted based on the parameters of the StepOnePlus™ 
Real-Time PCR System. 



Standard curves were prepared for each assay using 
known quantities of pCR™ 4 TOPO™ vector (Life Tech- 
nologies) containing the cloned human gene {ACTB, 120 
bp) and the 67 bp Leishmania infantum SSU rRNA frag- 
ment. The recombinant plasmids were serially diluted 
1:10 to create each standard curve. The quality param- 
eters of the standard curves, including PCR efficiency, 
linear dynamic range and correlation coefficient, were 
obtained by software analysis and were accurate and 
similar to those obtained in previous studies (Gomes 
et al. 2012, Marques et al. 2012). The parasite load was 
expressed by the Leishmania DNA load (relative copy 
number of the 67 bp SSU rRNA fragment) normalised 
against the reference gene ACTB, according to Over- 
bergh et al. (1999). ACTB copy numbers for the target 
samples were divided by the highest ACTB value ob- 
tained by the experiment, resulting in a correction factor 
used for normalisation. 

Statistical analysis - The Statistical Package for So- 
cial Sciences (SPSS®) 14.0 (IBM Corporation, New York, 
NY, USA) was used for data management and statistical 
analysis. Categorical variables were analysed by the chi- 
square test and, when appropriate, the Fisher exact test 
was used. For continuous variables, the Student t test 
and Wilcoxon test were utilised. The relative risk (RR) 
and the 95% confidence interval (CI) were calculated as 
a measure of effect of the dichotomous and categorical 
variables on the outcome. The variables with missing 
values above 50% were excluded from the analysis. A 
significance level of 5% was used. 

Ethics - This work was conducted in accordance with 
the Declaration of Helsinki and with Brazilian rules on 
human research. The study was approved by the Ethi- 
cal Committees of the Research Centre Rene Rachou/ 
Oswaldo Cruz Foundation and of FHEMIG. All parents 
or guardians signed an informed consent form prior to 
the study. 

RESULTS 

Forty-eight patients were included in the study. The 
average age was 3.3 years, ranging from three-nine years 
old. Thirty-seven children (77.1%) were under five years 
old and 25 (52.1%) were male. The median time between 
the onset of symptoms and hospital admission was 26 
days, ranging from four-180 days. 

In TO evaluation, the most frequent symptoms were 
as follows: fever (48 patients, 100%), pallor (42, 87.5%), 
malaise (39, 81.2%), abdominal distension (39, 81.2%), 
hyporexia (32, 66.7%), weight loss (27, 56.2%), vomit- 
ing (26, 54.2%) and abdominal pain (23, 47.9%). The 
most frequently observed signs on physical examination 
were the following: pallor (48, 100%), splenomegaly (47, 
97.9%) and hepatomegaly (46, 95.8%). 

The initial treatment was meglumine antimoniate in 
26 children (54.2%), amphotericin B deoxycholate in 15 
(31.2%), liposomal amphotericin B in two (4.2%) and a 
combination of a single dose of liposomal amphotericin 
followed by meglumine antimoniate in five (10.4%). In 
six patients (12.5%), the initial treatment was changed. 
In 37 patients, drug choice was based on assistant physi- 
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cian decision, according to Brazilian Ministry of Health 
recommendations. In other 11 patients, the initial drug 
was chosen based on a concomitant clinical trial that was 
taking place at the same hospital. Four patients who ini- 
tially received meglumine antimoniate finished the treat- 
ment with amphotericin B deoxycholate, due to an initial 
inadequate treatment response. In two other patients, 
amphotericin B deoxycholate was replaced by liposom- 
al amphotericin B, due to acute renal lesion. The mean 
treatment duration was 22 days. 



During hospitalisation, 16 patients (33.3%) met the 
criteria for poor outcome: 11 underwent blood transfu- 
sion, three underwent blood transfusion and were admit- 
ted to the ICU, one was admitted to the ICU and required 
mechanical ventilation and one underwent blood trans- 
fusion, was admitted to the ICU and required mechani- 
cal ventilation and vasoactive amine treatment. There 
was no death related to VL in the study period. 

The risk factors for poor outcome identified at TO were 
the following (Table): age < 12 months (RR 3.51; 95% CI 



TABLE 



Association between demographic, clinical and laboratory findings and poor outcome 
in children with visceral leishmaniasis hospitalised at Joao Paulo II Hospital, June 2010-June 2011 



Variables 


n 




Poor outcome 
(n) 

Yes No 


P 


Relative risk 
(95% CI) 


Age < 12 months 


48 


Yes 


0 


i 
i 


0.01" 


3.51 (1.89-6.52) 






No 


10 


31 






Symptoms > 60 days 


48 


Yes 


z 


"} 
J 


1.00" 


1.23 (0.39-3.90) 






No 


14 


z / 






Moderate or severe malnutrition 


48 


Yes 


1 
1 


Z 


1.00" 


1.00 (0.19-5.22) 






No 










Diarrhoea 


47 


Yes 




Q 


0.74" 


1.18 (0.49-2.82) 






No 


in 








Vomiting 


47 


Yes 


o 


1 7 

1 / 


0.66 


1.21 (0.51-2.86) 






No 


u 


I s 

1 J 






Bleeding (except cutaneous) 


48 


Yes 


1 


6 


0.40° 


0.39 (0.60-2.50) 






No 


15 


26 






Jaundice 


48 


Yes 


2 


0 


0.11" 


3.29 (2.12-5.09) 






No 


14 


32 






Oedema 


48 


Yes 


6 


5 


0.14" 


2.02 (0.95-4.30) 






No 


10 


27 






Tachydyspnoea 


48 


Yes 


3 


0 


0.03" 


3.46 (2.19-5.47) 






No 


13 


32 






Bacterial infection 


48 


Yes 


11 


9 


0.01 


3.08 (1.27-7.48) 






No 


5 


23 






Liver enlargement > 10 cm 


48 


Yes 


8 


4 


0.01" 


3.00 (1.44-6.23) 






No 


8 


28 






Spleen enlargement > 10 cm 


48 


Yes 


13 


19 


0.13 


2.17 (0.72-6.53) 






No 


3 


13 






Leukocytes < 1,500 cells/mm 3 


48 


Yes 


1 


2 


1.00" 


1.00 (0.19-5.22) 






No 


15 


30 






Neutrophils < 500 cells/mm 3 


48 


Yes 


3 


2 


0.32" 


1.98 (0.85-4.63) 






No 


13 


30 






Prothrombin activity < 70% 


41 


Yes 


10 


16 


0.74 


1.15 (0.49-2.74) 






No 


5 


10 






AST or ALT > 100 UI/L 


44 


Yes 


4 


10 


1.00" 


0.86 (0.32-2.26) 






No 


10 


20 






Total bilirubin > 1.0 mg/dL 


40 


Yes 


3 


2 


0.32" 


1.91 (0.80-4.54) 






No 


11 


24 






Albumin < 2.5 g/dL 


44 


Yes 


14 


8 


< 0.0001 


7.00 (1.80-27.24) 






No 


2 


20 






Creatinine clearance < 60 


43 


Yes 


1 


1 


1.00" 


1.58 (0.38-6.77) 


mL/min/1.73 m 2 




No 


13 


28 







a: Fisher's exact test; ALT: alanine aminotransferase; AST: aspartate aminotransferase; CI: confidence interval. 
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1.89-6.52), tachydyspnoea (RR 3.46; 95% CI 2.19-5.47), 
bacterial infection (RR 3.08; 95% CI 1.27-7.48), liver en- 
largement > 10 cm (RR 3.00; 95% CI 1.44-6.23) and se- 
rum albumin < 2.5 g / dL (RR 7.00; 95% CI 1.80-27.24). 

The qPCR was positive in all patients at TO. At Tl 
and T2, the parasitic loads were assessed in 46 and 43 pa- 
tients, respectively. Two children who were not assessed 
at Tl were assessed at T2. There was a wide variation in 
qPCR at each time point (Fig. 1). As expected, a statis- 
tically significant difference was observed between the 
number of SSU rRNA copies (in log 10 ) between T0-T1 (p 
< 0.0001) and T0-T2 (p < 0.0001), but no difference was 
observed between T1-T2 (p = 0.12). 

Among the 46 patients assessed at Tl, 35 (76.1%) had 
undetectable parasite DNA, in four (8.7%), SSU rRNA 
copies decreased and in the last seven (15.2%) the para- 
site load remained stable. Considering the 43 patients 
with a parasite load assessment at T2, 90.7% (39/43) had 
an undetectable parasite DNA; 29 of them were negative 
at Tl, nine were positive at Tl and became negative and 
one was not tested at Tl. Two children (4.7%) showed a 
decrease in the number of SSU rRNA copies (in log 10 ) 
compared to the previous evaluations, but the parasite 
DNA did not become undetectable (1 of them was not 
tested at Tl). Two other children (4.7%) had an undetect- 
able parasite DNA at Tl and a positive test at T2, with a 
number of copies similar to that at TO (Fig. 2). 

There was no statistical association between the num- 
ber of SSU rRNA copies (in log ]0 ) and the demographic, 
clinical and laboratory findings at TO. The mean parasite 
load was 2.94 SSU rRNA copies (in log 10 ) (95% CI 2.55- 
3.33) in the poor outcome group and 3.56 SSU rRNA 
copies (in log 10 ) (95% CI 3.07-4.06) in the other group. 
Despite a low p value, there was no statistical correlation 
between poor outcome and high parasite load at TO (p = 
0.08) (Fig. 3). 

DISCUSSION 

Age less than or equal to one year was related to in- 
creased lethality and poor response to pentavalent anti- 
mony use in other Brazilian studies (Santos et al. 2002, 
Rey et al. 2005, Costa et al. 2010). According to Braga et 
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Fig. 1: number of small-subunit ribosomal RNA (SSU rRNA) cop- 
ies (in log 10 ) before starting treatment (TO), between 10-15 days after 
initiation of treatment (Tl) and between 40-60 days after initiation of 
treatment (T2) in children with visceral leishmaniasis hospitalised at 
Joao Paulo II Hospital, June 2010-June 2011. 



al. (2013), presence of respiratory abnormalities at hos- 
pital admission was a risk factor for poor VL outcome in 
children. Tachydyspnoea was described as a risk factor 
for death by other authors (Costa et al. 2010, Sampaio 
et al. 2010) and its relation to disease severity can have 
multiple causes: interstitial pneumonitis secondary to 
VL, pneumonia or thoracic restriction due to increased 
sizes of the liver and of the spleen. Thus, its aetiology 
should be investigated carefully, including the use of 
X-ray images, to enable specific diagnosis and prompt 
treatment of the complication. 

The presence of bacterial infection in patients with 
VL is a recognised cause of death (Rey et al. 2005, Ol- 
iveira et al. 2006, Costa et al. 2010, Sampaio et al. 2010, 
Araujo et al. 2012) and was associated with poor out- 
come in this study. In FIJPII, febrile neutropaenic pa- 
tients (neutrophil < 500 cells/mm 3 ) with VL receive 
broad-spectrum antibiotics routinely, which could have 
positively influenced the disease course and may have 
contributed to the absence of deaths during the study 
period. 

Increased liver size was not associated with mortal- 
ity in previous studies. However, a Brazilian study sug- 
gested that the size of this organ is proportional to the 
duration of the illness (Pedrosa & da Rocha 2004). Thus, 
larger hepatomegaly may be associated with longer du- 
ration of the disease, which has already been associated 
with adverse VL outcome (Seaman et al. 1996, Collin et 
al. 2004). 

The reduction of albumin level, as observed in our 
study, has been associated with the severity of VL only 
in a study in Tunisia (Abdelmoula et al. 2003). Oedema 
and malnutrition, which are related to hypoalbuminae- 
mia and are associated with higher lethality (Collin et 
al. 2004, Rey et al. 2005), were not identified as risk fac- 
tors for poor outcome in this study. Therefore, once the 
change in albumin levels most likely precedes the clini- 
cal manifestations of malnutrition and oedema, monitor- 
ing of albumin levels may be relevant for early identifi- 
cation of potentially severe cases. 

All patients with VL showed a positive pre-treatment 
qPCR, indicating that it may be used as a diagnostic 
tool, in agreement with other studies in France and Italy 
(Mary et al. 2006, Antinori et al. 2009). This study also 
showed a good association between qPCR and treatment 
response: 76.1% (95% CI 60.9-87) showed an undetect- 
able parasite DNA during treatment and 90.7% (95% CI 
78.4-96.3) had an undetectable parasite DNA at the end 
of the treatment, which has also been described in stud- 
ies in France (Mary et al. 2004, 2006), Italy (Bossolasco 
et al. 2003), Tunisia (Aoun et al. 2009) and India (Sudar- 
shan et al. 2011). 

The four children who were positive by qPCR at T2 
were considered clinically cured and did not relapse dur- 
ing the period of the study, which suggests that parasitae- 
mia persistence or reemergence may not be associated 
with disease relapse, as described for immunosuppressed 
patients (Bossolasco et al. 2003, Mary et al. 2006). The 
positivity of qPCR at T2 could be explained by the detec- 
tion of free DNA instead of viable parasites (Fisa et al. 
2002) or by reinfection after return to endemic areas at T2. 
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The finding of two patients with negative qPCR at Tl and 
a positive test at T2 reinforces this hypothesis, but it can 
also be explained by a fluctuation in detectable DNA. 

There was no association between any of the symp- 
toms, signs, laboratory exams and the parasite load. This 
finding indicates that it is not possible to predict parasite 
load by any parameter of easier evaluation or lower cost 
than qPCR. 

This study was not able to show any association be- 
tween parasite load and clinical outcome. This finding 
differs from the results of previous studies in Sudan and 
Northeast Brazil (Seaman et al. 1996, Costa et al. 2010), 
which showed that a large number of parasites in the 
splenic aspirate and the BM aspirate, respectively, were 
related to a higher risk of death. These differences could 
be explained by the absence of association between the 



parasite counts in spleen, BM and peripheral blood and 
also by the different techniques used (optic microscopic 
vs. qPCR). The difference in study endpoints, death (Su- 
dan and Northeast Brazil) and poor outcome (current 
study) could also explain this difference. It is important 
to emphasise that no deaths were observed in this study. 
The wide variation in SSU rRNA copies (in log ]0 ), as 
reported by other authors in France (Mary et al. 2004, 
2006), Italy (Bossolasco et al. 2003), Tunisia (Aoun et al. 
2009) and India (Sudarshan et al. 2011) and the limited 
sample size may also have contributed to this finding. 

Although not statistically significant, there was a 
tendency of association between high parasite load and 
good outcome. This finding is in accordance with the 
high influence of individual immune status and other 
cofactors in disease manifestations (Murray et al. 2005). 
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Fig. 2: number of small-subunit ribosomal RNA (SSU rRNA) copies (in log |0 ) before starting treatment (TO), between 10-15 days after initiation 
of treatment (Tl) and between 40-60 days after initiation of treatment (T2) in 43 children with visceral leishmaniasis hospitalised at Joao Paulo 
II Hospital , June 2010-2011. A: 29 patients with an undetectable parasite DNA at Tl and T2; B: 10 patients with an undetectable parasite DNA 
only at T2; C: two patients with a decreased but detectable parasite DNA at T2; D: two patients with an undetectable parasite DNA at Tl and a 
detectable parasite DNA at T2. 
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Fig. 3: association between pre-treatment average number of small- 
subunit ribosomal RNA (SSU rRNA) copies (in log |0 ) and poor out- 
come among 48 children with visceral leishmaniasis hospitalised at 
Joao Paulo II Hospital, June 2010-June 2011. 



A higher parasite load might reinforce the specific im- 
mune response, which could enhance the host response 
and the disease resolution. But, considering the small 
sample size and large variation of parasite load, new data 
are needed to confirm these findings. 

The small number of patients is also a limitation of 
the present study, once it impaired a more extensive sta- 
tistical analysis, including multivariate analysis for poor 
outcome predictors. 

In this study, the definition of poor outcome may be a 
limitation, but it was necessary to redefine the endpoint 
once there are no deaths during the study period. Clinical 
outcome is influenced by the access to advanced medi- 
cal resources, including an ICU, mechanical ventilation 
and blood transfusions and by medical staff decisions 
and experience with the disease. The advantages of this 
study were the prospective design and the uniform as- 
sessment performed by a small number of doctors. 

In summary, age < 12 months old, tachydyspnoea, 
bacterial infection, liver size > 10 cm and serum albu- 
min < 2.5 g/dL were associated with poor outcome. Pa- 
tients with these pre-treatment findings should be care- 
fully monitored to offer proper assistance to those with 
a greater risk of complications. Pre-treatment parasite 
load in peripheral blood showed a good association with 
diagnosis and treatment response, but was not associated 
with poor outcome in paediatric patients. 
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